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Abstract: We present a new mechanism for baryogenesis, which links the baryon 
asymmetry of the universe to the dark matter density. The mechanism arises naturally 
in the Pentagon model of TeV scale physics. In that context, it forces a re-evaluation 
of some of the assumptions of the model, and we detail the changes that are required 
in order to fit observations. 



Contents 

1. Introduction 1 

2. Pseudo-goldstone dark matter with an asymmetry 2 

3. Spontaneous baryogenesis 4 

4. Spontaneous baryogenesis and dark matter in the Pentagon model 5 

5. Conclusions 8 

6. Acknowledgments 8 

7. Appendix 9 



1. Introduction 

With a few exceptions[l] most models of the early universe make no attempt to connect 
the observed baryon and dark matter densities. Dark matter is usually assumed to 
consist of neutralinos or axions, and there is no connection between the properties of 
these particles, and the baryon asymmetry. In this paper we will introduce a new class of 
models in which the solutions to these two problems are directly related. These models 
are motivated by the Pentagon model of TeV scale physics [2], where the mechanism 
we will discuss arises naturally. Indeed, it is forced on us if we insist on finding a dark 
matter candidate within the low energy model itself. 

The basic new structure in this class of models is an approximate symmetry, with 
current Jp B . This symmetry is explicitly broken and an asymmetry e for J PB is 
generated in the very early universe. The leading symmetry violating operator at low 
energy has dimension D and comes from physics at a scale The symmetry is 

also spontaneously broken at a scale /. Finally, there is a coupling in the low energy 
effective Lagrangian of the form 

— F B 

where is the ordinary baryon number current. 
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The asymmetry in Jpb then acts as an effective chemical potential for ordinary 
baryon number. Electroweak sphaleron processes will then generate a baryon asymme- 
try (and a compensating lepton asymmetry). This is a form of spontaneous baryogen- 
esis, a mechanism invented by Cohen and Kaplan[3]. 

At a lower energy scale, the effects of spontaneous and explicit violation of Jpb 
set in. These convert the energy stored initially in the asymmetry into a condensate of 
pseudo-Goldstone bosons, which can be the dark matter. We can use the parameters 
of the model to fit both the baryon asymmetry, and dark matter density, in a variety 
of ways. However, in the version of the model that arises from the Pentagon, all 
parameters but e and are determined in terms of other experimental quantities. We 
find that, given the assumption in [2] that M b is at least as large as the neutrino seesaw 
scale (5 x 10 14 GeV)the baryon asymmetry is too large, or the dark matter density too 
small. We describe modifications of the model which remove this problem in section 4. 

2. Pseudo-goldstone dark matter with an asymmetry 

The idea of pseudo-Nambu-Goldstone bosons as dark matter is familiar from axion 
models. The general idea is that, after inflation, the pseudo Goldstone field has a fixed 
value and small velocity. Its motion is friction dominated until the Hubble parameter 
falls to the axion mass, at which time it begins to oscillate and behaves like cold dark 
matter. For appropriate values of the mass of the field and its decay constant, a PNGB 
can reproduce the observed dark matter density. 

Here we propose an alternative initial condition for PNGB dark matter. The 
symmetry associated with the PNGB is explicitly broken. We assume that this breaking 
is larger during inflation than it is at low energies. A large value for the inflaton field, 
typical in slow roll inflation models, can enhance operators that are highly irrelevant 
at low energies. We assume that this produces an asymmetry e at the end of inflation, 
and that there is not much entropy production between the time this asymmetry is 
generated, and the present. Note that the Affleck-Dine mechanism can easily generate 
large values of e[5]. 

We note parenthetically that in string theory most PNGB's arise as duals of fun- 
damental anti-symmetric tensor gauge fields. From this point of view, the assumption 
of an asymmetry is equivalent to assuming an expectation value for the magnetic field 
strength and has been recently studied by Ibanez[4]. Our work will be primarily 
devoted to PNGBs of accidental symmetries of low energy gauge fields, but for com- 
pleteness we will also study what happens to a QCD axion with asymmetric initial 
conditions. 
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The non-zero value of e and our assumption of isentropic cosmological expansion 



allows us to write an equation for the time dependence of the density J° 



PB 



J°pb = e dT 3 , (2.1) 
where g is the effective number of massless degrees of freedom. This equation is valid 
from the moment of asymmetry generation, until some time at which low energy Jp B 
violation becomes important, or a large amount of entropy is dumped into the universe. 
Note that it is valid both above and below the scale of spontaneous breaking of the 
approximate symmetry. Below that scale, which we denote by /, we have 

4b = fdob, (2.2) 

where b is the PNGB field. The primary focus of this paper is on models in which / is 
dynamically generated by strong gauge dynamics. 

Assume that the explicit breaking of the symmetry after inflation comes primarily 
from an operator of dimension D, and has its origin in physics which decouples at the 
scale Mb. Below the scale / this will give rise to an effective potential for b 

%iV(blf). (2.3) 



M b 



For a QCD axion we would have D = 4 and = Aqcd- The cross over temperature, 
at which this potential begins to affect the evolution of the b field is given by 

{eg{T*Yf ph D L 



M fe D - 4 ' 



(2.4) 



or 

A D/2 

eg(T*) 3 = V2fM^^- (2.5) 
At this temperature, the ratio of energy densities of dark matter to radiation is: 

-Pi 



P61 to(T*) 3 ) 2 



(2.6) 



r* f 2 g(T*y ■ 

Below this temperature, the b field begins to oscillate, and its energy scales like that of 
cold dark matter. The dark matter to radiation ratio at any lower temperature is 

PjL (eg(T*)*)* T* 

An asymmetric PNGB can be a good dark matter candidate if this ratio becomes one 
at matter radiation equality, T eq ~ leV. Using the equation for T*, this requires 

Mr 2 = -^t a l 2 - (2-8) 

J 1 eg 
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3. Spontaneous baryogenesis 



Suppose in addition that there is a coupling 

Y 2 JpbB„ (3.1) 

between the current of the PNGB and the ordinary baryon number current. Assume 
for the moment that T* < T sh) where T sh is the scale below which the electroweak 
baryon violating process shuts off, with T sh ~ 100 GeV. Then, in the regime where 
electroweak baryon violation is in equilibrium, this coupling has the form 

eg^B = fiB , (3.2) 

which is a time dependent chemical potential for baryon number. The relative rate of 
change of the chemical potential is the expansion rate of the universe, much slower than 
sphaleron processes. Thus, the combination of jx and the sphaleron process put us in the 
regime of equilibrium spontaneous baryogenesis, as defined by Cohen and Kaplan[3]. 
Other baryon violating processes could be hypothesized, at a variety of energy scales. 
Their contribution to the total baryon asymmetry would add to the one we compute 
here. Without fine tuning, it is unlikely that the contributions of other processes could 
cancel the electroweak sphaleron induced asymmetry. It is straightforward to compute 
the induced baryon asymmetry (see Appendix). We will also assume that any chemical 
potential for lepton number (recall that B — L is conserved by sphaleron processes) is 
much smaller than that for baryon number, and that B — L asymmetries generated 
in the early universe are small compared to the baryon asymmetry generated by our 
mechanism. We find 

The last equality is the constraint from the observed baryon asymmetry. 

Note that if T* is substantially larger than T s h then no asymmetry is generated 
because the chemical potential is turned off before baryon number violating processes 
go out of equilibrium. The asymmetry will equilibrate to whatever value of the B — L 
asymmetry was generated by early universe processes like leptogenesis. 

Let us plug this value of e into the formula we got by insisting that the dark matter 
density comes out right, then we get 

Mr 2 = 1.2 x 1(T 9 -^— A? /2 . (3.4) 

/ 1 sh 1 eq 
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The condition that T* < T s h translates into the inequality 



Af < gM b D ~ A T! h T eq: 



(3.5) 
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(3.6) 



while the plausible constraint that e < 10 3 i 



is 



1.2(— ) 2 < 10 



(3.7) 



sh 



For a QCD axion we expect A > /, where a strictly greater than sign is used because the 
coupling between the axion field and the baryon number current violates CP and might 
be suppressed by more than dimensional analysis. The conventional lower bounds on 
the axion decay constant, from red giants and supernovae, then rule out this kind of 
asymmetric axion scenario. Note also that such a scenario would have required a value 
of e which is probably too large to be generated by the Affleck-Dine mechanism. 

Our analysis was based on the assumption that the pseudo-Goldstone boson was 
the correct description of physics at the scale where the primordial Jp B asymmetry is 
wiped out by processes which violate this symmetry. The temperature T* where this 
occurs must thus be smaller than /. Scenarios where this inequality is not satisfied are 
more complicated to analyze. Some of them could give rise to acceptable cosmologies. 
However, both the QCD axion models, and the Pentagon model, satisfy T* < f so we 
will not attempt to analyze this possibility any further. 

4. Spontaneous baryogenesis and dark matter in the Pentagon 



In the Pentagon model, all of the previous ingredients are present, and most of the 
parameters are related. There is a spontaneously broken accidental symmetry, penta- 
baryon number. It has f — y TeV and A 5 = x TeV, where x and y are of order 
1 — 10 and most probably y > x Using various different estimates of standard model 
superpartner masses that have appeared in the literature, we find x running between 
1.5 and 7. There is also a current current coupling to baryon number with A5 = CK3A, 
where the strong coupling is evaluated at the TeV scale and is ~ .1. 

The lowest dimension penta-baryon number violating operators, which preserve 
gauge invariance, SUSY, and the fundamental discrete R symmetry of the model are 



model 




(4.1) 
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and 

f d 2 9 SP b . (4.2) 



These have dimension D = 7 1 . In [2] TB also imposed the natural constraint that all 
allowed irrelevant operators were suppressed by powers of Mjj ~ 10 15 GeV, which is the 
scale that appears in the neutrino seesaw term. Let us let the scale in the dimension 7 
operator be a free parameter, M b . 

The requirement that we get the right baryon asymmetry is 

£B = l e( nt )2 - (4 - 3) 

or 

e = 1.2x 2 x 1(T 5 . (4.4) 



The equation determining T* is 



{eg{T*Yf _ Aj 
2/ 2 ~M b 3 ' 



(4.5) 



or 



A 7/2 

e 9 {T*f = yftf-jLp. (4.6) 

According to WMAP, the temperature of matter radiation equality is about 1 eV. Thus, 
the condition that the penton is dark matter is 

Pb {eg{T*ff T* 

p 7 Pg{T*Y T eq [J) 
If we insert the values for e and T eq into this equation, we get the constraint 

M b = (|^) 2/ V x 10 8 GeV. (4.8) 

The prefactor is a number slightly less than 1, so this ranges between about 3 x 10 8 — > 
3x 10 10 GeV as x ranges over the values allowed by the various estimates of superpartner 
masses. We also note that over the whole range of x, the inequality T* < T sh is satisfied, 
as long as y/x is not too large. 

The low scale M b raises the specter of unacceptably fast proton decay. There is a 
dimension 6 operator of the schematic form 

^2 f d20 Q* LS ' 



1 In [2] TB forgot the R symmetry constraint, and used the D = 6 operators without the singlet S. 
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which is invariant under all of the symmetries of the Pentagon model. With M = Mj, 
this would lead to disaster. Thus, if we want to use the penton to generate the dark 
matter in the universe, we must construct a theory at the scale M b which explains the 
absence of operators of dimension 6, which could contribute to proton decay. 

The alternative is to abandon the penton field b as the origin of dark matter, which 
is to say that e is presumed to be small. We could retain a non-zero value of e as the 
mechanism for baryogenesis, though this would be no more attractive than a host of 
other options. Dark matter would have to come from somewhere outside the Pentagon. 
A possible candidate is an axion dual to an antisymmetric tensor field. This could also 
solve the strong CP problem and be compatible with estimates from string theory[7] . 

Yet a third alternative, is to raise the values of A 5 and m IS s, probably abandoning 
the hypothetical connection of the model to Cosmological SUSY Breaking. A value 
of A 5 ~ 3 x 10 5 GeV and M b = 10 15 GeV seems to be compatible with both the 
observed dark matter density and baryon asymmetry. The value of e is a bit less 
than 1 and T* is well below the sphaleron mass. However, there are a number of 
phenomenological problems with this suggestion. The rough estimates for sparticle 
masses and the electroweak scale in the Pentagon model are 

m i/2 ~ ~3 9S J^ miSS) ( 9 ) 




/50 Oil , t -,r,\ 

H u ~ 240GeV ~ yA 5 . (4.11) 

j^m, ss . (4.12) 

miss is the mass term which creates a meta-stable SUSY violating vacuum in the 
Pentagon model. The value of A runs between 1 and 8, while that of B is between 
1 and 4. gs is a Yukawa coupling between a next-to-minimal SSM singlet, and the 
Pentagon fields. Other squark and slepton masses are larger than that of the right 
handed selectron by factors of 

The factor of 1/6 in the equation for the electroweak scale represents a pious hope 
that the original Pentagon model with A 5 ~ 1.5 TeV does not suffer from a little 
hierarchy problem. That is, the dimensional analysis estimate of the electroweak scale 
has 1/6 replaced by 1, and we have to hope that dynamical calculations in the strongly 
coupled Pentagon model provide the factor of 6. 
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If we postulate A 5 ~ 3 x 10 5 GeV, dynamical suppression is no longer plausible. 
We can get the correct electroweak scale by choosing g s small, about 5 x 1CT 3 , but this 
implies small gaugino masses. The ratio between squark and wino masses is 

ma A CHq . 
VTLyj IB a 2 

Using the experimental lower bound on the wino mass we get squark masses that 
are a few times 10 TeV . The model then has a hierarchy problem, and radiative 
corrections to the Higgs mass are substantially larger than the values indicated by 
precision electroweak fits. 

Even a value M b = 10 15 GeV is not enough to protect us from proton decay. 
The unified values of gauge couplings are quite large in the Pentagon model, so even 
dimension six proton decay operators must be suppressed by 10 16 GeV or so. We 
conclude that raising the scale A 5 does not seem to be a promising avenue for making 
a working model of penton dark matter and baryogenesis. 



5. Conclusions 

The Pentagon model suggests a novel form of spontaneous baryogenesis, which can tie 
together the dark matter density and the baryon asymmetry. While the general idea 
works quite well, it does not work in the Pentagon model, unless we contemplate a 
scale of 10 8 — > 10 10 GeV for the leading irrelevant operator that violates penta-baryon 
number. It remains to be seen whether one can invent a high energy extension of the 
model to generate this operator without generating standard model operators already 
ruled out by experiment. 

If we stick with 10 15 GeV as the scale for all irrelevant corrections to the Lagrangian, 
then we must abandon the penton theory of dark matter. The penton is a light PNGB 
which might be produced in the laboratory, but its cosmological abundance is negligible. 
It could still participate in the generation of the baryon asymmetry we observe. 

The most plausible candidate for dark matter compatible with the Pentagon would 
then be a QCD axion. This would also solve the strong CP problem. 
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7. Appendix 

In this appendix we go through the calculation of the baryon number expectation, in 
the presence of a chemical potential in thermal equilibrium. We begin by ignoring 
leptons, which would play a role in determining the equilibrium baryon number gener- 
ated through processes such as the electroweak sphaleron where the quantity B — L is 
conserved. After computing (B) first by ignoring this constraint, we will then modify 
our calculation to include it. This is mostly for pedagogical purposes, but the naive 
calculation turns out to give the right order of magnitude as long as the lepton chemical 
potential is of equal or smaller magnitude than the baryon chemical potential. 

Statistically, the equilibrium value of B can be calculated from the partition func- 
tion: 

The n being used here is dimensionless. In terms of the usual chemical potential jl, 
which has dimensions of energy, /i — ^ — {1(3. 

The partition function for the baryons including the chemical potential is 

1 — Tr{e~^ H+li $ dAxq \ qi ) 

where 



i,k 

1 X 

d 4 xq t q= - Y^WUk ~ b{b k ] 



k 

This can be re-expressed as 

Z = Tr[exp(- - |)afa?>) exp(- + 

k % , k 

= Tr[HeM-(P\k\ - f )af a») exp(-(/?|ifc| + ^)bfhf)] 

i,k 

The trace gives a sum over all possible states of the system; for fermions, each state 
can only have occupation number equal to zero or one, thus 



Z = J](l + e -WI-f))(l + e-WI + t)) 



k.i 
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We are interested in the logarithm of Z, which allows us to express the logarithm of 
the product as a sum of logarithms 

lnZ = 6x3x2x J^K 1 + e~ (/3|fe| +t) + ln(l + e - (m ~^)] 
k 

The factor of 6 x 3 x 2 comes from the sum over % which includes six flavors of quarks, 
three colors, and two spin states. Taking k to the continuum limit we have 



and taking the /i derivative gives 

d , „ „„ r d 3 k 



(B) = — lnZ = 36V / 



(27r)3 L 3(e ^+^)) 3(^(*-&)) J 



Using spherical coordinates, the angular integral just gives a factor of 4n, which leaves 
an integral over the magnitude k: 

6V 



7T 2 



where 



h = / dk- 



-I 

-L 



e«*-6> + 1 



OO K 



e/*(*+&) + i 



By substituting x = f3{k — ^), I\ becomes: 



3/3 V 

dx- ; 



/3 3 J_i£ e* + l /3 3 7 e^ + l /9 3 7_i£ + 1 



and by substituting x = (5{k + ^), J 2 becomes: 



/ 2 = 



? 3 7 t e* + l /3 3 7 e* + l (3 3 J 



x-jy 

f3 3 Je ~~ + 1 /3 3 y + 1 f3 3 J ~~ + 1 

The two can then be recombined into: 

i r r . MS) f° , ( x + 1) 2 r° , 



h ~ h ~J 3 



r dx ^ + r dx^t+fdx^m 

Jo e x + l J_e e x + l Je e x + 1 J 



3 3 
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h-h 



1 



)(-)+ / dxC— ' 



A/i n 2 



+ 



/ e x + 1 e x + 
Miraculously, the exponentials simply add up to 1, leaving an integral of a polyno- 
mial: 



1 



rc 2 fi 
~9~ 



+ 



dx{x 



Q / 



-[-(-) + -(-) 3 ] 

[3 3[ 3 V 3 V J 



(B) 



l 



3/? 3 ^ 27n 2 ' 

The above discussion, however, was naive in that we have ignored the role of 
leptons, ie this calculation does not conserve B — L. Therefore, we must modify the 
original partition function by including a chemical potential for leptons similar to that 
for baryons, and impose this constraint by including a Kronecker delta function: 



1 — T r [ e -P H e VB f d 4 xq}q ie iJ, L f d 4 xl}hfi B ^ 



where now 



i,k 

with the Baryon and Lepton operators expressed as 

J d A xq ] q = i Y^la{a k - b[b k ], 

d A xlH = y^[<4c fe - d\d k }. 



The delta function can be represented in integral form as 

5bl = ±- f dae^ B -V 
lix 



(We can choose any range covering a full period of 2-7T for the integration limits, but as 
we will find later, the range from — n to n turns out to be particularly convenient.) 

B-L = E(5« ) -i« i, -4 i M i, + « ) ) 
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It follows that 



Z = Tr — 




As before, the trace gives a sum over all possible states of the system (zero or one 
for fermions), thus: 



Unlike the case for baryons alone, we cannot exploit the fact that we wish to find 
the logarithm of Z to convert this product into a sum of logarithms because of the a 
integral. However, we can write 



Again we have taken k to the continuum limit. As before the baryons obtain a 
factor of 36 from the sum over % (6 flavors, 3 colors, 2 spin states), the leptons similarly 
obtain a factor of 3 generations x 3 spin states (two for each charged lepton and one 
each neutrino) = 9. 

The integral in the exponent is the same integral we had in the purely baryonic 
case, but with the real parameter /x/3 replaced by a complex parameter z, where 
z = (//£ + ia)/3 for the baryonic terms, and z = hl — ia for the leptonic terms. The 
integrand has poles along the imaginary axis of the z-plane at Im{z} = n + 2nn, n 
integer, and they occur when k = Re{z}. The analytic form of the integral we found 
for real z above, 



is valid as long as Imj^;} is restricted to the range (— For Imjz} outside of 
this range, the function simply repeats with period 2n, which happens because the 




k,i k,i 




+9[ln(l + e (-/»l*l+w--*0) + ln(l + e (-/^|-^+*a))] 
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integrand was periodic in a. The closed form expression for the full periodic function 
can be written by replacing z with hie 2 in the polynomial. As long as we choose the 
range of integration over a to be from — n to n, the original polynomial form above can 
be used. The partition function is then: 



The constant term is independent of both //# and /il- Since we are only interested in 
expectation values which do not depend on the normalization of the partition function, 
we can ignore this term. 

The a integral can be performed using the saddle point approximation. Allowing 
a to be complex now, there is a maximum in Re{/} at Re{a} = 0. The saddle point 
occurs at the point on the imaginary axis where there is a minimum in Re{/} along the 
imaginary direction in the a-plane. In the text, we will find that /il,b must be small 
in order to fit the observed baryon asymmetry. Therefore, we will work to lowest order 
in the chemical potentials, and avoid writing the lengthy exact expression. 




where 



7T 2 1 7T 2 1 

f(a,HB,ftL) = 36[— ((^B + ia)/3) 2 + — ((/i£ + ia)/3) 4 ] +9[— (^L-ia) 2 + — (fi L -ia) 4 } 
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To first order in hl and /j,b, the saddle point is at: 



4i 9i 

^saddle — Y3^ B ~~ 13^ L 

The value of / at caddie, to second order in \ib and \il is: 

/(laddie, Hb,Hl) = 9.11 / u b /!l + 4.55/4 + 4.55^ 

Z can be approximated well by evaluating the integrand at a sa ddie- In the saddle 
point approximation, this would usually need to be divided by the square root of the 
second derivative of / with respect to a, but this just leads to an additive term in 
InZ which is independant of the volume of spacetime. As V — > oo, this is negligible 
compared to the term proportionate to V. 



111 Z ~ 2^/^ asaddle ' ^ = 0M ]p( fX BtiL + + + 0{n BL ) 

Therefore, the expectation values of baryon and lepton number are, to first order 
in \ib and 

._. ._. 91nZ 91nZ V , . 

{B) = {L) = l^ = l^* i W {fiB + fiL) 

This calculation gives us the baryon asymmetry written in Equation 3.3, when we plug 
in the effective chemical potential in Equation 3.2. 
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